A Simple Varying-alpha Cosmology 
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We investigate the cosmological consequences of a simple theory in which the electric charge e is 
allowed to vary. The theory is locally gauge and Lorentz invariant, and satisfies general covariance. 
We find that in this theory the fine structure 'constant', a, remains almost constant in the radiation 
era, undergoes small increase in the matter era, but approaches a constant value when the universe 
starts accelerating because of the presence of a positive cosmological constant. This model satisfies 
geonuclear, nucleosynthesis, and CMB constraints on time-variation in a, while fitting simultane- 
ously the observed accelerating universe and the recent high-redshift evidence for small a variations 
in quasar spectra. The model also places specific restrictions on the nature of the dark matter. 
Further tests, involving stellar spectra and the Eotvos experiment, are proposed. 



There is renewed interest in cosmological theories 
where the traditional constants of Nature can vary in 
space and time. Theoretical interest is motivated by 
string and M-theories in which the true constants exist in 
more than 3 + 1 dimensions and the effective (3+ 1)— di- 
mensional constants can display cosmological variations 
in time and space. The most obscrvationally sensitive 
'constant' is the electromagnetic fine structure constant, 
a = e 2 /he, and recent observations motivate the for- 
mulation of varying-a theories. The new observational 
many-multiplet technique of Webb et al, 13], pf, exploits 
the extra sensitivity gained by studying relativistic tran- 
sitions to different ground states using absorption lines 
in quasar (QSO) spectra at medium redshift. It has pro- 
vided the first evidence that the fine structure constant 
might change with cosmological time [[j] || . The trend of 
these results is that the value of a was lower in the past, 
with Aa/a = -0.72±0.18x lCT 5 for z « 0.5-3.5. Other 
investigations have claimed preferred non-zero values of 
Act < to best fit the cosmic microwave background 
(CMB) and Big Bang Nucleosynthesis (BBN) data at 
z fa 10 3 and z fa 10 10 respectively [j|||. 

A varying a may be closely related to a varying speed 
of light (VSL), [||-J8), which is of interest because it ap- 
pears to solve the cosmological problems resolved by in- 
flation, together with some other problems M . VSL the- 
ories are usually associated with breaking of Lorentz in- 
variance (but see |ll[). A less radical approach is to 
attribute varying a to changes in the fundamental elec- 
tron charge, e. Bekenstein p2f gives an example of a 
consistent varying-e theory which preserves local gauge 
and Lorentz invariance, and is generally covariant. This 
is a dilaton theory with coupling to the electromagnetic 
"F 2 " part of the Lagrangian, but not to the other gauge 
fields. The approach adopted is only a matter of con- 
venience in the choice of units: they are observationally 
indistinguishable |7],|lC],|l3| . 

Another remarkable set of recent observations is of 
Type la supernovae in distant galaxies. These data 
have extended the Hubble diagram to redshifts, z > 1 
p3. They imply an accelerated expansion of the uni- 



verse. When combined with CMB data, the super- 
novae observations favour a flat universe with approxi- 
mate matter density, fl m fa 0.3 and vacuum energy den- 
sity, Q\ fa 0.7. Studies have attempted to determine 
whether quintessential scalar fields could explain both 
cosmological dark matter and the recent acceleration of 
the universe, |15-19|. 

We will not here attempt to explain the acceleration 
of the universe. Instead, we show that by applying a 
generalisation of Bekenstein's varying-e theory in a cos- 
mological setting including the cosmological constant, A, 
we are able to explain the magnitude and sense of the ob- 
served change in a. The main assumption is that the cold 
dark matter has magnetic fields dominating their electric 
fields. The magnetostatic energy then drives changes in 
a in the matter dominated epoch, but as the Universe 
starts to accelerate these changes become friction domi- 
nated and come to a halt. This gives a decelerated rate 
of change in a, just as the universe starts to accelerate, 
in accord with both data sets. The only energy scale 
we introduce is of the order of Planck scale, which also 
makes our model attractive. This model may be seen as 
a more conservative alternative to |2(j,|35|], where a VSL 
scenario was proposed which could explain the observed 
acceleration of the universe and variations in a, as well 
as their remarkable coincidence in redshift space. 

Bekenstein's original theory takes c and h to be con- 
stants and attributes variations in a to changes in e, or 
the permittivity of free space. This is done by letting 
e take on the value of a real scalar field which varies in 
space and time cq —> e = eoe(a; M ), where e is a dimension- 
less scalar field and eo is a constant denoting the present 
value of e. This means some well established assump- 
tions, like charge conservation, must give way pi| . Still, 
the principles of local gauge invariance and causality are 
maintained, as is the scale invariance of the e field. 

Since e is the electromagnetic coupling, the e field cou- 
ples to the gauge field as eA^ in the Lagrangian and the 
gauge transformation which leaves the action invariant is 
eA^ — ► eA^ + x,fi, rather than the usual A^ — > A^ + x, M - 
The gauge-invariant electromagnetic field tensor is then 



Ffiv = (( e ^),^ — (e-4/j),i/) A) which reduces to the usual 
form when e is constant. The electromagnetic Lagrangian 
is still C em — —F^Ffiv/A and the dynamics of the e field 
are controlled by the kinetic term C e — — ^u>(e t fj,e'^)/e 2 , 
as in pa] (we use a metric with signature -+++). Here 
the coupling constant u) can be written as jt, where I is 
the characteristic length scale of the theory, introduced 
for dimensional reasons. This constant length scale gives 
the scale down to which the electric field around a point 
charge is accurately Coulombic. The corresponding en- 
ergy scale, hc/l, has to lie above a few tens of MeV to 
avoid conflict with experiment. 

Consider Bekenstein's theory in the cosmological set- 
ting suggested by the recent supernovae results. To sim- 
plify calculations, we invoke a transformation introduced 
in ref. p2 ] . By defining an auxiliary gauge potential a^ = 
eA^, and field tensor f^ = eF^ v — d^a^ — d v a^, the co- 
variant derivative takes the usual form, D^ = d^ + ieoa^. 
The dependence on e in the Lagrangian then occurs only 
in the kinetic term for e and in the F 2 = f 2 /e 2 term. To 



simplify further we change variable: e 
total action becomes 
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where C+ = -%d^d^ and C em = -\f^f» u . The 
gravitational Lagrangian is the usual C g — 1 „ G R 1 with 
R the curvature scalar. Our theory generalises Beken- 
stein's approach by including the effects of the varying e 
(or ip) field on the gravitational dynamics of the expand- 
ing universe. The scalar field ip plays a similar role to the 
dilaton in the low-energy limit of string and M-theories, 
with the important difference that it couples only to elec- 
tromagnetic energy. Since the dilaton field couples to all 
the matter then the strong and electroweak charges, as 
well as particle masses, can also vary with x^ . These sim- 
ilarities highlight the deep connections between effective 
fundamental theories in higher dimensions and varying- 
constant theories, p3[ . 

To obtain the cosmological equations we vary the ac- 
tion with respect to the metric to give the generalised 
Einstein equations 
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and vary with respect to the ip field to give the equations 
of motion for the field that carries the a variations: 
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(3) 



It is clear that C em vanishes for a sea of pure radiation 
since then C em — (E 2 — B 2 )/2 = 0. This suggests a 
negligible change in e in the radiation epoch, a fact con- 
firmed by our numerical calculations. The only signifi- 
cant contribution to a variation in ip comes from nearly 
pure electrostatic or magnetostatic energy. 

In the matter epoch changes in e will occur. In or- 
der to make quantitative predictions we need to know 
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FIG. 1. Cosmological evolution from radiation domination 
through matter domination and into lambda domination with 
coupling C m /dj = —0.02%. The upper graph plots a as a 
function of the redshift z. The lower graph shows the energy 

densities of radiation ( ), dust ( ), cosmological constant 

( ) and the scalar field (combined) as fractions of the total 

energy density. The scalar field energy is subdominant at 
all times, a increases in the matter era, but approaches a 
constant after A takes over the expansion. 



how non-rclativistic matter contributes to the RHS of 
Eqn. (||). This is parametrised by the ratio £ = C em /p, 
where p is the energy density, and for baryonic matter 
C em « E 2 /2. For protons and neutrons ( p and C, n can 
be estimated from the electromagnetic corrections to the 
nucleon mass, 0.63 MeV and —0.13 MeV, respectively 
p0| . This correction contains the E 2 /2 contribution (al- 
ways positive), but also terms of the form j^a^ (where 
jp is the quarks' current) and so cannot be used directly. 
Hence we take a guiding value ( p rs ( n ~ I0~ 4 . Further- 
more the cosmological value of (, (denoted C, m ) has to be 
weighted by the fraction of matter that is non-baryonic, 
a point ignored in the literature ll2U24]. Hence, Cm de- 
pends strongly on the nature of the dark matter, and it 
could be that (cdm w — 1 (e.g. superconducting cos- 
mic strings, for which C em ~ —B 2 /2), or |Ccda/ | *C 1 
(neutrinos). BBN predicts an approximate value for the 
baryon density of fls ~ 0.03 with a Hubble parameter 
of ho ~ 0.6 , implying flcDM ~ 0.3. Hence, depending 
on the nature of the dark matter, £ m can be positive or 
negative and have a modulus between and s=s 1. 

Assuming a spatially- flat, homogeneous and isotropic 
Friedmann metric with expansion scale factor a(t) we 
obtain the Friedmann equation (G = c = 1) 
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where the cosmological vacuum energy p\ is a constant 
given by A/(87r), and p^ — jip 2 - For the scalar field we 
get 




FIG. 2. The data points are the QSO results for the chang- 
ing a(z) reported in refs. [pJ-0]. The solid line is the theoret- 
ical prediction for a(z) in our model with £m/a> = —0.02%. 

The top ( ) and bottom ( ) lines correspond to choices 

Cm/u) = —0.01% and ^ m /u — —0.03% respectively 
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where H = a/a. The conservation equations for the non- 
interacting radiation and matter densities, p r and p m 
respectively, are: 
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0, with 



This last relation can be written as p r + AHp r 
p r = p r e~ 2 ^ ex a~ 4 . 

Eqns. (0-0), govern the Friedmann universe with time- 
varying a = exp(2i/))eg/?ic. They depend on the choice of 
the parameter £ m /w, which we take to be negative. We 
evolve these equations numerically from early radiation- 
domination, through the matter era and into vacuum 
domination by pa- Fig. n] shows the evolution of a with 
redshift in this model, for (, m /w = —0.02%. We note that 
ip and a remain almost constant during early radiation 
domination where baryonic species become relativistic. 
In the matter epoch, ip and a increase slightly towards 
lower redshifts, but tend to constant values again once 
the universe starts accelerating, and A dominates - this 
is due to the friction term Hip in Eq. (||). This A effect 
reduces variations in a during the last expansion time of 
our universe where the local geonuclear effects of varying 
a are strongly constrained by observations, |2rj , |27|, while 
allowing the cosmological variations observed by [hi El at 
redshifts, z k, 0.5 — 3.5, where the effects of A on the 
expansion progressively diminish. In Figure we plot 
the predicted change in a for —C, m JLO — 0.01, 0.02, 0.03%, 
and the binned QSO data from refs. jj]-S|. Given the un- 
certainties in C, m discussed above, it is possible to fit the 



data with w = 0(1), so that the theory's length scale is 
of the order of the Planck length. 

In view of the special a(z) variation produced by the 
cosmic acceleration there is agreement with all labora- 
tory, geological and astrophysical constraints on varying- 
a deriving from the last expansion time (cf. [pi|-p7|). No- 
tice also that the supernovae luminosity data are fitted 
by our model, since ip affects the cosmological expansion 
very little, and its direct effect upon the luminosities of 
astrophysical objects is negligible. Hence, our Hubble 
diagram is precisely the same as that of a universe with 
constant a and O m w 0.3 and fl\ ~ 0.7. Our model also 
meets constraints from BBN, since it occurs deep in the 
radiation epoch, z ~ 10 9 — 10 10 , when a is predicted to 
be only 0.007% lower than today. The standard BBN 
scenario can withstand variations in a of the order of 1% 
without contradicting observations (see W and references 
therein). The value of a at CMB decoupling, z « 1000 is 
only ~ 0.005% lower than today, compatible with CMB 
observations, QH, which place an upper bound of a few 
percent. However, the variations we predict are close 
enough to these limits to hold out the possibility of obser- 
vational test in the future by more detailed calculations 
of the effects on BBN and the CMB, and more precise 
data. Low-redshift observations of molecular and atomic 
transitions ||28| can provide important information about 
the value of a close to the redshift where acceleration 
commences, z\ ~ 0.7, if the chemical isotopic evolution 
uncertainty can be reduced p9[ . 

Spatial variations of a are likely to be significant J52| , 
and our model makes firm predictions on how a varies 
near massive objects. Linearising eq. (J3J) and following 
the calculation of [B3|, one finds that relative variations 
in a are proportional to the local gravitational potential: 
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where M is the mass of the compact object, r is its ra- 
dius, and £, is its value of £. When £ m and £ s have dif- 
ferent signs, for a cosmologically increasing a, we predict 
that a should decrease on approach to a massive object. 
If |Cm| ~ Cs, on extragalactic scales the CMB tempera- 
ture anisotropy AT/T ~ GM/r would lead us to expect 
large-scale spatial gradients of order Aa/a ~ 10 -9 . More 
locally, one would need an object not larger than some 
tens of Schwarzschild radii for the effect on a(r) to be 
observable with current technology. However with im- 
proved technology, one might find less demanding candi- 
dates. An independent low-z test of the effects seen by 
plpl could ultimately be provided by the detection of a 
Aa 7^ effect from the fine structure of stellar spectral 
lines. The exact relation between the change in a with 
redshift and in space (near massive objects) is model de- 
pendent p3| , but eq.(@) provides the exact prediction 
for the simple varying-a theory considered in this pa- 
per. Notice that these variations do not conflict with the 
Pound- Rebka-Snider experiments ]34[|. 



Spatial gradients in a lead to an extra force acting 
upon matter coupling to V y i a the f^ m term. Since 
Cp ¥" (n this force acts differently on matter with different 
composition leading to violations of the weak equivalence 
principle pflmt] . These are parameterized by the Eotvos 
parameter, which in our theory is 
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where E denotes the Earth and 1 and 2 two different 
test bodies. If £ n rs C P ~ ICp — Cn| the first factor is 
O(10 -5 ) for typical substances used in experiments. The 
third factor, C, m /to, is of the order of — 10 -4 . Hence for 
I Cm I = 0(0.1) — O(l) we have consistency with the cur- 
rent experimental bound, rj < 10~ 13 [p4j| . We note that 
the next generation of Eotvos experiments should be able 
to detect the variations in a predicted by this theory, but 
firmer predictions require better theoretical calculations 
of £ for neutrons, protons, nuclei and atoms (the uncer- 
tainties of which were discussed above) . 

In summary, we have shown how a cosmological gener- 
alisation of Bckenstein's theory of a varying e can natu- 
rally explain the reported variations in the fine structure 
constant whilst satisfying all other observational bounds. 
The onset of A domination is shown to be closely related 
to the cosmic epoch when significant changes in a cease 
to occur. Our numerical results show that with a natu- 
ral coupling, and using observational constraints on the 
nature of the cold dark matter, a changes significantly 
only in the matter dominated epoch. At the onset of 
A domination, the expansion accelerates and a rapidly 
approaches a constant. This model also places specific 
restrictions on the nature of the dark matter. 

Our model complies with geonuclear constraints, like 
Oklo, but is consistent with the non-zero variations in 
a{z) inferred from observations of quasar absorption lines 
[Q-g[ at z « 0.5 — 3.5. It is attractive because apart from 
the (observed) cosmological constant value, the only free 
parameter introduced is an energy scale of the order of 
the Planck scale. It is also consistent with CMB and 
BBN observational constraints. Further tests for this 
model will be possible using stellar spectra and the next 
generation of Eotvos experiments. 
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